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Abstract 

We study diagonal multipoint Pade approximants to functions of the 
form 

Fi.) = f "^^ + Riz), 



where i? is a rational function and A is a complex measure with compact 
(-H regular support included in E, whose argument has bounded variation on 

the support. Assuming that interpolation sets are such that their normal- 
ized counting measures converge sufficiently fast in the weak-star sense 
Q to some conjugate-symmetric distribution a, we show that the counting 

' ' measures of poles of the approximants converge to a, the balayage of a 

. onto the support of A, in the weak* sense, that the approximants them- 

K^ selves converge in capacity to F outside the support of A, and that the 

VQ poles of R attract at least as many poles of the approximants as their 

^^ multiplicity and not much more. 

CN AMS Classification (MSC2000): primary 41A20, 41A30, 42C05; secondary 

^ ^ 30D50, 30D55, 30E10, 31A15. 
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nomials, non-Hermitian orthogonality. 

1 Introduction 

This paper is concerned with the asymptotic behavior of diagonal multipoint 
Pade approximants to functions of the form 

n.)-/^+i?(.), (1.1) 

where i? is a rational function holomorphic at infinity and A is a complex measure 
compactly and regularly supported on the real line. 

Diagonal multipoint Pade approximants are rational interpolants of type 
(n, n) where, for each n, a set of 2n-|- 1 interpolation points has been prescribed, 
one of which is infinity. Moreover, we assume that the interpolation points con- 
verge sufficiently fast to a conjugate-symmetric limit distribution whose support 
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is disjoint from both the poles of R and the convex hull of supp(A), the support 



of A (see (|2.8|)). 

To put our results into perspective, let us begin with an account of the ex- 
isting literature. When A is a positive measure and _R = (in this case F is 
referred to as a Markov function), the study of diagonal Pade approximants 
to F at infinity goes back to A. A. Markov who showed (see [53]) that they 
converge uniformly to F on compact subsets of C \ /, where / is the convex 
hull of supp(A). Later this work was extended to multipoint Pade approximants 
with conjugate-symmetric interpolation schemes by A. A. Gonchar and G. Lopez 
Lagomasino in |16j . A cornerstone of the theory is the close relationship between 
Pade approximants to Markov functions and orthogonal polynomials, since the 
denominator of the n-th diagonal approximant is the n-th orthogonal polyno- 
mial in L^{dX) (resp. L^{dX/p), where p is a polynomial vanishing at finite 
interpolation points). For further references and sharp error rates, we refer the 
reader to the monographs |32l ES] ■ 

Another generalization of Markov's result was obtained by A. A. Gonchar 
on adding polar singularities, i.e. on making R ^ 0. He proved in [T5] that 
Pade approximants still converge to F locally uniformly in C \ (5" U /), where 
S" is the set of poles of R, provided that A is a positive measure with singular 
part supported on a set of logarithmic capacity zero. Subsequently, it was 
shown by E. A. Rakhmanov in [57] that weaker assumptions on A can spoil the 
convergence, but at the same time that if the coefficients of R are real, then the 
locally uniform convergence holds for any positive A. Although it is not a concern 
to us here, let us mention that one may also relax the assumption that supp(A) 
is compact. In particular, Pade and multipoint Pade approximants to Cauchy 
transforms of positive measures supported in [0, oo] (such functions are said to 
be of Stieltjes type) were investigated by G. Lopez Lagomasino in [THI^D]. Let 
us also stress that polynomials satisfying certain Sobolev-type orthogonality 
exhibit an asymptotic behavior quite similar to that of the denominators of 
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The case of a complex measure was taken up by G. Baxter in [TI] and 
by J. Nuttall and S. R. Singh in [53], who established strong asymptotics of 
non-Hermitian orthogonal polynomials on a segment for measures that are ab- 
solutely continuous with respect to the (logarithmic) equilibrium distribution 
of that segment, and whose density satisfy appropriate conditions expressing, 
in one way or another, that it is smoothly invertible. These results entail that 
the Pade approximants to F converge uniformly to the latter on compact sub- 
sets of C \ / when R = and dX/dfii meets these conditions (here fij indi- 
cates the equilibrium distribution on /). For instance Baxter's condition is that 
log dX/dfij, when extended periodically, has an absolutely summable Fourier 
series. When d\{t)/dt is holomorphic and nonvanishing on a neighborhood of /, 
still stronger asymptotics, which apply to multipoint Pade approximants as well, 
were recently obtained by A. I. Aptekarev in [4J (see also 5 ), using the matrix 
Riemann-Hilbert approach pioneered by P. Deift and X. Zhou (see e.g. [12]). 
Eventhough it is not directly related to the present work we mention for com- 
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pleteness another approach to analyzing the asymptotics of Pade approximants 
based on three term recurrence relations |10j . 

Meanwhile H. Stahl opened up new perspectives in his pathbreaking papers 
[3TJ|31j, where he studied diagonal Pade approximants to (branches of) multiple- 
valued functions that can be continued analytically without restriction except 
over a set of capacity zero (typical examples are functions with poles and branch- 
points). By essentially representing the "main" singular part of the function as 
a Cauchy integral over a system of cuts of minimal capacity, and through a deep 
analysis of zeros of non-Hermitian orthogonal polynomials on such systems of 
cuts, he established the asymptotic distribution of poles and subsequently the 
convergence in capacity of the Pade approximants on the complement of the 
cuts. In [17] this construction was generalized to certain carefully chosen mul- 
tipoint Pade approximants by A. A. Gonchar and E. A. Rakhmanov, who, in 
particular, used it to illustrate the sharpness of O. G. Parfenov's theorem (for- 
merly Gonchar's conjecture) on the rate of approximation by rational functions 
over compact subsets of the domain of holomorphy, see |25]. Of course the true 
power of this method lies with the fact that it allows one to deal with measures 
supported on more general systems of arcs than a segment, which is beyond the 
scope of the present paper. However, since a segment is the simplest example 
of an arc of minimal logarithmic capacity connecting two points, the results 



we just mentioned apply in particular to functions of the form (1.1 1, where A 
is a complex measure supported on a segment which is absolutely continuous 
there with continuous density that does not vanish outside a set of capacity 
zero. By different, operator-theoretic methods, combined with a well-known 
theorem of E. A. Rakhmanov on ratio asymptotics (see [28j), A. Magnus fur- 
ther showed that the diagonal Pade approximants to F converge uniformly on 
compact subsets of C \ / when R = and dX/dt is non-zero almost everywhere 
with continuous argument |22j . The existence of a uniformly convergent subse- 



quence of diagonal Pade approximants to ( 1.1 1 with non-trivial R was shown in 
|34] whenever supp(A) is a disjoint union of analytic arcs in "general position" 
of minimal capacity and dX/dt is sufficiently smooth and non- vanishing. More- 
over, when supp(A) is a union of several intervals and the density of the measure 
is real analytic, the behavior of the zeros that do not approach supp(A) nor the 
poles of R can be described by the generalized Dubrovin system of non-linear 
differential equations [55] . 

In contrast with previous work, the present approach allows the complex 
measure A to vanish on a large subset of /. Specifically, we require that the 
total variation measure |A| has compact regular support and that it is not too 
thin, say, larger than a power of the radius on relative balls of the support 
(see the definition of the class BVT in Section l2l. In particular, this entails 
that supp(A) could be a thick Cantor set, or else the closure of a union of 
infinitely many intervals; such cases could not be handled by previously known 
methods. Although fairly general, these conditions could be further weakened, 
for instance down to the A-criterion introduced by H. Stahl and V. Totik in 
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Y\ However, our most stringent assumption bears on the argument of A, as 
we require the Radon-Nikodym derivative (iA/(i|A| to be of bounded variation 
on supp(A). This assumption, introduced in [HI [6], unlocks many difficulties 
and will lead us to the weak convergence of the poles and to the convergence in 
capacity on C \ (S" U supp(A)) of multipoint Pade approximants to functions of 



the form (1.1 1. Moreover we shall prove that each pole of R attracts at least as 
many poles of the approximants as its multiplicity, and not much more. In fact, 
our hypotheses give rise to an explicit upper bound on the number of poles of 
the approximants that may lie outside a given neighborhood of the singular set 
of F. Hence, on each compact subset K oi C\{S' U supp(A)), every sequence 
of approximants contains a subsequence that converges uniformly to F locally 
uniformly on K \ E, where E consists of boundedly many (unknown) points. 
When supp(A) is a finite union of intervals, results of this type were obtained 
under stronger assumptions in [25 for classical Pade approximants. 

Finally, we would like to mention that the presented approach can also be 
carried out for AAK-type meromorphic approximants. Although their defini- 
tion is rather simple, deriving functional decomposition for them is not trivial 
(cf. [T] and [S]) and the latter gives rise to more comphcated orthogonality re- 
lations than those satisfied by the denominators of Pade approximants. Thus, 
we consider meromorphic approximants separately in [S]. 

2 Pade Approximation 

We start by describing the class of measures that we allow in (fOJ and placing 
restrictions on the points with respect to which we shall define Pade approxi- 
mants. 

Let A be a complex Borel measure whose support S :— supp(A) C M is 
compact and consists of infinitely many points. Denote by | A| the total variation 
measure. Clearly A is absolutely continuous with respect to |A|, and we shall 
assume that its Radon-Nikodym derivative (which is of unit modulus |A|-a.e.) 
is of bounded variation. In other words, A is of the form 

dA(i) = e'"^(*)d|A|(i), (2.2) 

for some real-valued argument function (p such thalrl 



Viifi, S) :- sup <^ V \ipix,) - (^(x,_i)| } < oo, (2.3) 




where the supremum is taken over all finite sequences xq < xi < . . . < xn in S 
as N ranges over N. 

For convenience, we extend the definition of (p to the whole of M as follows. 
Let / := [a, b] be the convex hull of S. It is easy to see that if we interpolate 



^This depends on the corresponding generalization of the results in 6_ to be found in |18| . 
as yet unpublished. 

■^Note that e'"^ has bounded variation if and only if ip can be chosen of bounded variation. 
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ip linearly in each component of I \ S and if we set ip{x) := linit_+a^ tes v(^) for 



X < a and ip{x) :— limt^f,^ tes 'fii't) ior x > b (the limits exist by (2.3 1), the 
variation of if will remain the same. In other words, we may arrange things so 
that the extension of ip, still denoted by ip, satisfies 

Vip,S) = V{cp,R)^:V{ip). 



Among all complex Borel measures of type (2.2l-(2.3l, we shall consider 
only a subclass BVT defined as follows. We say that a complex measure A, 
compactly supported on M, belongs to the class BVT if it has an argument of 
bounded variation and if moreover 

(1) supp(A) is a regular set; 

(2) there exist positive constants c and L such that, for any x G supp(A) and 
S G (0,1), the total variation of X satisfies \X\{[x — d,x + 6]) > cS . 

In what follows we consider only functions of the form 

F{z):^J^ + R.{z), (2.4) 

with A G BVT and Rg a rational function of type (s — 1, s) assumed to be in 
irreducible form. Hereafter we shall denote by 

Qs{z)^l[{z-vr^^-''> (2.5) 

the denominator of Rg, where S' is the set of poles of Rg and m{ri) stands for 
the multiplicity of 77 G 5'. Thus, i^ is a meromorphic function in C \ S with 
poles at each point of S" and therefore it is holomorphic in C\ S, where 

S:=SUS'. 

Note that F does not reduce to a rational function since S consists of infinitely 
many points (cf. j7] Sec. 5.1] for a detailed argument). 

Diagonal multipoint Pade approximants to F are rational functions of type 
(n, n) that interpolate _F at a prescribed system of points. More precisely, pick 
n G N and let A^ = {Ci.td ■ • • j C2n,n} be a set of 2n interpolation points, where 
the Q,n G C \ 5' need not be distinct nor finite. With such an A^ we form the 
monic polynomial 

V2n{z)= n (^-O.n) (2.6) 

(note that V2n retains only the interpolation points at finite distance thus it 
needs not have exact degree 2n). 



Given F of type (2.4-) and An as above, the diagonal multipoint Fade approxi- 
mant to F associated with An is the unique rational function n„ = Pn/Qn where 
the polynomials pn and g„ satisfy: 
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(i) degpn < n, degg„ < n, and q„ ^ 0; 
(ii) {qn{z)F{z) -p„(z)) /v2n{z) is analytic inC\ S; 

(iii) {qn{z)F{z) - p,,{z)) /V2n{z) = O {1/Z^+^) aS Z ^ ^ . 

A multipoint Pade approximant always exists since the conditions for p„ and 
Qn amount to solving a system of 2n+l homogeneous linear equations with 2n+2 
unknown coefficients, no solution of which can be such that qn = (we may 
thus assume that qn is monic); note that (iii) entails at least one interpolation 
condition at infinity and therefore n„ is, in fact, of type {n — l,n). 

If we let now A := {An}ni£N be an interpolation scheme, i.e. a sequence 
indexed by n S N of sets An as above, we get a corresponding sequence {n„}„gN 
of diagonal Pade approximants whose asymptotic behavior can be studied when 
n gets large. Namely, we shall be interested in three types of questions: 

(a) What is the asymptotic distribution of the poles of Pade approximants to 
F? 

(b) Do some of these poles converge to the polar singularities of F? 

(c) What can be said about the convergence of such approximants to F? 

To be able to provide answers to these questions, we need to place some 
constraints on interpolation schemes. An interpolation scheme A is said to be 

admissible if 

(1) /C(A), the set of the limit points of A, is disjoint from S' U /; 

(2) the counting measures of the points in A^ converge in the weak* topology 
to some Borel measure, say a, having finite logarithmic energy; 



(3) the argument functions of polynomials U2n, associated to A via (2.6), have 
uniformly bounded derivatives on /. 

In other words, we call an interpolation scheme admissible if the interpolation 
points stay away from the poles of R^ and the convex hull of the support of A, 
if there exists a Borel measure a = cr(A) supported on /C(A) such that 

_. 2n 

n ^-^ ^' 

and if the norms ||(w2n/|'y2n|)'||/ are uniformly bounded with n, where || • \\k 
stands for the supremum norm on a set K. We call a the asymptotic distribution 
of A. Note that /C(A) is not necessarily compact. If it is not compact, the 
finiteness of the logarithmic energy of a is understood as follows. Since /C(A) is 
closed and does not intersect S, there exists zq € C\UkAk such that zg ^ /C(A). 
Pick such a zq and set Mz^iz) := l/(z — zq). Then, all Mza{Ak) are contained in 
some compact set and their counting measures converge weak* to a^ such that 
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a'^{B) := a{M-^{B)) for any Borel set -B C C. We say that A is admissible if 
a'' has finite logarithmic energy. Obviously, this definition does not depend on 
a particular choice of zq. Further, as a consequence of (3), there exists a finite 
constant Va satisfying 

V{&Tg{v2n),I)<VA for any n G N. (2.7) 

Notice that (3) is satisfied if, for example, all A„ in A are conjugate-symmetric. 
More generally, it can be readily verified that (3) amounts to 

uniformly on /, which is exactly what we meant in the introduction when saying 
that the counting measures of interpolation points should converge sufficiently 
fast. 

The four theorems stated below constitute the main results of the paper. For 
the notions of potential that we use (logarithmic and Green potentials, balayage, 
equilibrium distributions, capacity and convergence in capacity) the reader may 
want to consult the appendix. 



Theorem 2.1. Let F be given by (24)-(2.5) with A G BVT and let {Iin}nefi 
be a sequence of diagonal multipoint Fade approximants to F that corresponds 
to an admissible interpolation scheme A with asymptotic distribution a. Then 
the counting measures of the poles of n„ converge in the weak* sense to a, the 
balayage of a onto S. 

We note that the limit distribution of poles of n„ can also be interpreted as 
the weighted equilibrium distribution on S in the presence of the external field 
-U" (cf. [30, Ch. I]). 

Recall (cf. |30|, pg. 118]) that S^o is simply fis, the logarithmic equilibrium 
distribution on S. Therefore for classical Fade approximants (when each V2n = 
1, i.e. when all the interpolation points are at infinity), the above theorem 
reduces to the following result. 



Corollary 2.2. Let F be given by (24)-(2.5) with A e BVT and let {n„}„gN be 



the sequence of Fade approximants to F at infinity. Then the counting measures 
of the poles of Yin converge to fis in the weak* sense. 

The previous theorem gave one answer to question (a). Our next result 
addresses question (c) by stating that the approximants behave rather nicely 
toward the approximated function, namely they converge in capacity to F on 

C\S'. 



Theorem 2.3. Let F, A, and {n„}„gN be as in Theorem 2.1 Then 



|(F-n„)(z)|V2"-Pexp{-C/cV(^)} (2.9) 

on compact subsets ofC\S, where U^xg is the Green potential of a relative to 



cap 



C\ S and — > denotes convergence in capacity. 
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Finally, we approach question (6). In order to provide an answer to this 
question, we need some notation. For any point z G C define the lower and 
upper characteristic m(z) , rn{z) G Z+ as 

fn( z) :— iid Wi( z , U) , fn(z,U) :— lini max#{5„nC/|, 

U AT-^oo n>N 



and 



m(z) := inf m(z, C/), m(z,U) := lim Toin ^{SnfMJ}, 

U N^oo n>N 



respectively, where the infimum is taken over all open sets containing z and 
Sn is the set of poles of n„, counting multiplicities. Clearly, rn{z) < m(z). 



m(z) = +CX) if 2 g 5 by Theorem 2.1 and m(z) = if and only if z is not a 
limit point of poles of n„. Further, let /,„ := {[o-j, ^j]}i'Li be any finite system 
of intervals covering S. Also, let Arg(^) e (— tt, tt] be the principal branch of 
the argument, where we set Arg(O) = tt. With this definition, Arg(-) becomes a 
left continuous function on M. Now, for any interval [flj, bj] in /„ we define the 
angle in which this interval is seen at f € C by 

Angle(C, [aj,b,]) := |Arg(aj - C) - Arg(&, - 01- 

Finally, we define additively this angle for the whole system, i.e. the angle in 
which /,„ is seen at ^ is defined bjjj 

m 

0(O:=^Angle(e,K-,6,]). (2.10) 

Note that < 9{^) < n and 6'(^) = tt if and only if ^ e /„. 

The forthcoming theorem implies that each pole of F attracts at least as 
many poles of Fade approximants as its multiplicity and not much more. 



Theorem 2.4. Let F, A, and {Iln}neN be as in Theorem 2.1 and 6{-) be the 
angle function for a system of m intervals covering S. Then 

m{ri) > m{ri), rj € S' , (2-11) 

and 

J2 {rn{v)~Hv)){^-0{v))<V, (2.12) 

neS'\S 

with 

V ■.= V{<fi) + VA + {m + 2s'-l)Tr + 2 ^ m{ri)9{i]), (2.13) 

rieS'\S 



where Va was defined in {2.1) and s' is the number of poles of R on S counting 
multiplicities. 



^The notation does not reflect the dependency on the system of intervals, but the latter 
will always be made clear. 
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The basis of our approach lies in analyzing the asymptotic zero distribution 
of certain non-Hermitian orthogonal polynomials. It is easy to understand why. 
Indeed, let T be any closed Jordan curve that separates S and /C(A) and contains 
S in the bounded component of its complement, say D. Since 

iqnF~Pn){z)/v2niz)=0{l/z''+') aS Z -^ OO 

and the left-hand side is analytic in C \ S*, the Cauchy formula yields 

r dz 

/ z^q,Xz)F{z) ^-0, j = 0,...,i 



, n-1, zeD. 



rn,k 



It V2n{z) 

Clearly, by writting Rg as 

771(77)— 1 

))GS' fe=0 ^ " 

we see that the last equations are equivalent to 



/"-("'"("H-E s 



fc! V Vonii) 



= (2.14) 



for all Pn-i G 7^n-i by the definition of F, the Fubini-Tonelli's theorem, and the 
residue formula. So, upon taking P„_i to be a multiple of Qs, these relations 
yield for n > s 

t''Qs{t)qn{t)^^ = 0, k = 0,...,n-s-l. (2.15) 

Hence the denominators of the multipoint Fade approximants to F are poly- 
nomials satisfying non-Hermitian orthogonality relations with varying complex 
measures dX/v2n- 

The following theorem describes the zero distribution of the polynomials g„ 



satisfying (2.15). Let us stress that, in general, such polynomials need not be 
unique up to a multiplicative constant nor have exact degree n. In the theorem 
below, it is understood that g„ is any sequence of such polynomials and that 
their counting measures are normalized by 1/n so that they may no longer be 
probability measures. This is of no importance since the defect n — deg(g„) is 
uniformly bounded as will be shown later. 

Theorem 2.5. Let {qn}neN be a sequence of polynomials of degree at most 



n satisfying weighted orthogonality relations (2.15), where {v2n}neN is the se- 
quence of monic polynomials associated via (2.6) to some admissible interpo- 
lation scheme A with asymptotic distribution a and where X € BVT. Then 
the counting measures Vn of the zeros of qn{z) = JK-^ ~ 0,")? namely i/„ := 
(l/n) ^ (5j. „, converge in the weak* sense to a, the balayage of a onto S = 
supp(A). 
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By virtue of the results in the PhD thesis of R. Kiistner [TB] , a generalization 
of the previous theorem can be proved when the measure A, instead of belonging 
to BVT, has an argument of bounded variation and satisfies the so-called A- 
criterion introduced in |33l Sec. 4.2]: 

f(. o 1- Log(l//i[i -r,i + r]) n 

cap {teS: hmsup — ^^ '^'^ , , \ < +oo } = cap(S^). 

Vl '-^0 Log(l/r) jy 

However, this assumption would make the exposition heavier and we leave it to 
the interested reader to carry out the details. 

3 Proofs 

We start by stating several auxiliary results that are crucial for the proof of 
Theorem 1231 



Lemma (fSJ Lem. 3.2]) Let v he a positive measure which has infinitely 
many points in its support and assume the latter is covered by finitely many 
disjoint intervals: supp(i^) C U"L]^[aj, fej]. Let further ^j be a function of bounded 

variation on supp(i/). // the polynomial ui{z) — IliLi(^ ~ Cj); di < I, satisfies 

f t^ui{t)e'''^*Uv{t) = 0, fc = 0, . . . , ? - 1, 

then 

di m 

Y,{n ~ e{^,)) + {l~ di)n < Y, Vi^, K-, &,]) + (m - l)7r, 



where 0{-) is the angle function defined in (2.10) for a system of intervals 
ur=i[a,,6,]. 

As a consequence of this lemma, we get the following. 

Lemma 3.1. Let qn{z) = 11 7=1 (-^^ ~ ^j.n) ^6 ^^^ n-th orthogonal polynomial in 



the sense of (2.15), where X S BVT and the polynomials U2n are associated to 
an admissible interpolation scheme A. Then 

£('^-^(0.n)) + ("-c?n)^< V(</7) + Fa+^ m(77)%) + (m + s-lK, (3.16) 
i=i ves' 



where Va was defined in [2.1) and 9{-) is the angle function defined in (2.10) 
for a system of intervals Im '■— U"li[aj, bj] that covers S with I — [ai, bm] being 
the convex hull of S. 

Proof: Denote by i^n{t) an argument function for e^^''*'^Qs{t)qn{t) /v2n{t) on 
/, say 

ipnit) = v{t) ~ arg(z;2„(i)) + ^ m(77)Arg(i - ry) + ^ Arg(i - ^,,„). 

JjGS' i=l 
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It is easy to see that tpn is of bounded variation. Further, set I = n — s, 

>s{t)qn{t) 



ip = ipnj and, diy(t) 



V2n{t) 



d\xm. 



Then it follows from orthogonality relations (2.151 that 
t^e''l'^'Uv{t) = 0, fc = 0, , 



, n — s — 1. 
Thus, the previous lemma, applied with 7i; = 1, implies that 

m 

y^ V{ipn, [aj.bj]) >{n-s-m + l)7r. 
So, we are left to show that 

i = l rjeS' i=l 

By the definition of '0n, we have 

771 777 771 

777 

+ mZ ™('/)1^(Arg(- - ?7), K-, &j]) 

m dn 

+ EE^(Arg(--C^n),K-,6,]). 

3=1 7=1 

The assertion of the lemma now follows from the fact that, by monotonicity, 
y(Arg(.-0,K6])-Angle(e,[a,6]). 



Finally, we state the last two technical observation. 

Lemma 3.2. With the previous notation the following statements hold true 

(a) Let ^ he a real function of hounded variation on an interval [a, h] and 
Q a polynomial. Then there exists a polynomial T ^ and a constant 
(3 € (0,7r/32) such that 



Arg U'''^''^Q{x)T{x)) < 7r/2 - 2/3 
for all X G [a, h] such that T{x)Q{x) ^ 0. 



(3.17) 
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(b) Assume that the polynomials V2n o,re associated to an admissible inter- 
polation scheme. Then for every e > there exists an integer I and a 
polynomial Ti^n of degree at most I satisfying: 



V2nix) 

Tl^n{x) 



|'y2n(a:)| 



< e, X G I, 



for all n large enough. In particular, the argument of Ti^n/v2n l^es in the 
interval (— 2e, 2e) for such n. 

Proof: (a) When Q = I this is exactly the statement of Lemma 3.4 in [3] 
and since ip{x) + Arg((3(x)) is still a real function of bounded variation on /, 



( |3.17| follows. 

(5) This claim follows from Jackson's theorem [13, Thm. 6.2] since the deriva- 
tives of V2n/\v2n\ ^^^ Uniformly bounded on /. ■ 
Note that Lemma |3.1| applied with to — 1, implies that the defect n — dn is 
bounded above independently of n. 

Corollary 3.3. Let U be a neighborhood of S . Then there exists a constant 
kjj £ N such that each g„ has at most ku zeros outside of U for n large enough. 

Proof: Since U is open, its intersection with ( — 1, 1) is a countable union of 
intervals. By compactness, a finite number of them will cover 5*, say U"L]^ (oj ,bj). 
Apply Lemma [3?T] to the closure of these intervals intersected with / and observe 
that any zero of (?„ which lies outside of U will contribute to the left-hand side 



of (3.161 by more than some positive fixed constant which depends only on U. 



Since the right-hand side of (3.161 does not depend on n and is finite we can 
have only finitely many such zeros. ■ 

Proof of Theorem \2.3( Observe that we may suppose A is contained in a compact 
set. Indeed, if this is not the case, we can pick a real number xq ^ ^[^(A) U S" U / 
and consider the analytic automorphism of C given by Mxo{z) := l/{z — xq), 
with inverse M-1{t) = a;o + 1/t. If we put A^ := M^„{A.a), then A' = {A^} 
is an admissible interpolation scheme having asymptotic distribution a\ with 
a^{B) = a{AI^^^{B)) for any Borel set B C C Moreover, the choice of Xq yields 
that IC{A^) is compact. Now, if we let 

inir) - tV(M-/(t)), 

L,(r) = r^Qs{M-^{T)), 

PLs-iir) - r— ip„_,_i(M-i(r)), 

vLir) - T'-V2n{M-!{r)), 

then £„ is a polynomial of degree n with zeros at M^^^^jn), J = 1, • ■ • , dn, and a 
zero at the origin with multiplicity n — dn. In addition, W2„ is a polynomial with 
a zero at each point of A^^, counting multiplicity. Thus, up to a multiplicative 



constant, Ujn i^ the polynomial associated with A^^ via (2.6 1. Analogously, L 



is a polynomial of degree s with a zero of multiplicity m(ri) at Mx^iri), rj € S' , 



L.Baratchart and M.Yattselev 13 



and Pn-s-1 is ^'^ arbitrary polyno mial of degree at most n — s — 1. Making the 
substitution t — M~^^{t) in (2.151, we get 



where dX'^{T) = rdX [M~^{t)) is a complex measure with compact support 
Mx:g{S) C M, having an argument of bounded variation and total variation 
measure |A'*| G BVT. Note that r is bounded away from zero on supp(A''), since 
S is compact and therefore bounded away from infinity. Now, since Lemma |3.1| 
implies that n — dn is uniformly bounded above, the asymptotic distribution of 
the counting measures of zeros of in is the same as the asymptotic distribution 
of the images of the counting measures of zeros of qn under the map M^o ■ As 
the counting measures of the points in A^^ converge weak" to cr" , it is enough to 
show that counting measures of zeros of £„ converge to ct*^ since the balayage 
i s pres erved under Mxo(e.g. because harmonic functions are, cf. equation 
(A. 61 1 in the appendixj^J Hence we assume in the rest of the proof that A is 



contained in a compact set, say Kq, which is disjoint from S by the definition 
of admissibility. 

Now, let r be a closed Jordan arc such that the bounded component of 
C \ r, say D, contains S while the unbounded component contains Kq. Then 

Qn = qn,l ■ qn.2, whcrC 

QnA^) = n i^-^j.n) and g„^2(^) = JJ (z - ^j>). (3.18) 



Corollary 3.3 assures that degrees of polynomials qn.2 are uniformly bounded 
with respect to n, therefore the asymptotic distribution of the zeros of g„^i 
coincides with that of g„. Denote by ly^.i the zero counting measure of g„ i 
normalized with 1/n. Since all Vn,! are supported on a fixed compact set, Helly's 
selection theorem and Corollary |3 . 3| yield the existence of a subsequence Ni such 
that I'n.i -^ V for n G Ni and some Borel probability measure v supported on 
S\ remember the defect n — deg(g„^i) is uniformly bounded which is why :/ is a 
probability measure in spite of the normalization of g„_i with 1/n. 

Next, we observe it is enough to show that the logarithmic potential oiv~a 
is constant q.e. on S. Indeed, since supp(cr) is disjoint from S and subsequently 
t/~'^ is harmonic on 5, U" is bounded q.e. on S under this assumption. Hence, 
by lower semi-continuity of potentials, C/" is bounded everywhere on 5* and 
therefore v has finite energy. The latter is sufficient for v to be C-absolutely 
continuxyuaj Moreover, we also get in this case that V"" is constant q.e. on 



5* by (A. 60 1 and, of course, u is also C-absolutely continuous. Thus, v ~a hy 



the second unicity theorem [30!, Thm. II. 4. 6]. 

*Here we somewhat abuse the notation and use the symbol '^ to denote the balayage onto 
l^xoiS), while in the rest of the text it always stands for the balayage onto S. 

*A Borel measure /x is called C-absolutely continuous if fJ,{E) = for any Borel polar set. 
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Now suppose that W^ '^ is a constant q.e. on S. Then there exist nonpolar 
Borel subsets of S, say E^ and E^ , and two constants d and r > such that 



U''-''{x)>d + T, xeE+, V^^ix) <d-2T, xeE^ 
Then we claim that there exists j/o G supp(i') such that 

[/"--"(yo) > d. 

Indeed, otherwise we would have that 
U''{x)<U"{x) + d, a;esupp(j/). 



(3.19) 



(3.20) 



Then the principle of domination [301, Thm. II. 3. 2] would yield that (3.20) is 



true for all z G C, but this would contradict the existence of E^. 

Since /C(A) is contained in the complement of D, the sequence of potentials 
{C^'^"}nGNi converges to V^ locally uniformly in D. This implies that for any 
given sequence of points {j/n} C D such that i/n ^ yo 8£ n ^ cx3, n e Ni, we 
have 



hm U''-{y^)^U''{yo). 



(3.21) 



On the other hand, by applying the principle of descent [30, Thm. 1.6.8] for the 
above sequence {?/„}, we obtain 



liminf C/^^-Hyn) > C^'(2/o). 



Combining (3.191, (3.21), and (3.22) we get 



liminf C/'^-^~""(2/„)>C/''-"(z/o)>d. 



(3.22) 



(3.23) 



Since {j/n} was an arbitrary sequence in D converging to yo, we deduce from 
( 3.23 1 that there exists p > such that, for any y E [yo — 2p, j/o + 2p] and n G Ni 
large enough, the following inequality holds 



Clearly 



U" 



'(y) 



1 
2n 



log 



V2n{y) 



liiiy) 



(3.24) 



(3.25) 



and therefore inequality (3.24) can be rewritten as 



grli(y) 

V2n{y) 



<e-'"^ ye[yo-2p,yo + 2p]. 
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for all n S Ni large enough. We also remark that the same bound holds if 
{<ln,i} is replaced by a sequence of monic polynomials, say {un}, of respective 
degrees n + o{n), whose counting measures normalized by 1/n have asymptotic 
distribution v. Moreover, in this case 



qn,iiy)uniy) 



V2n{y) 



< e 



-2nd 



(3.26) 



for any y & [yo — 2/3, yo + 2p] and all n G Ni large enough. 

In another connection, since U'^~'^{x) < d — 2t on E-, applying the lower 
envelope theorem [30l Thm. 1.6.9] we get 

liminf [/''"•I"'"" (a;) = V^^^ix) < d - 2t, for q.e. x(^E^. (3.27) 

n — *oo, nGNi 

Let Z be a finite system of points from /, to be specified later, and denote for 
simplicity 

K{z) = ql^i{z)/v2n{z). 

Then by [H [3] there exists So C S such that 5*0 is regular, cap (-E_ n Sq) > 
and dist(Z, S'o) > 0, where dist(Z, S'o) := min^g^ dist(z, 5'o). Thus, there exists 
X € E^ n So such that 

|6„(a;)|>e-2"('i-), n e N2 C Ni, 



by (3.25) and (3.27). Let Xn be a point where |6„| attains its maximum on Sq, 



M„ 



||?>«||so = l&n(x„)l>e-2"''~"^. 



(3.28) 



Since f2n has no zeros in D, the function log |6„| is subharmonic there. Thus, 
the two-constant theorem [251 Thm. 4.3.7] on 13 \ 6*0 yields 



log|6„(z)| <log(M„) LUD\So{z,So) + 2nlog (^r^^^^ 



{l~UJD\Soiz,So)), 



z € D, where w^jx Sp is the harmonic measure on D\Sq, d(D) := max{diam(£)), 1}, 
and diam(Z3) := xnaxx_yi=D \x ~ y\- Then we get from (3.28) that 

l-uio\Soiz,So) f r\iD\ \ 2"(l-"D\So(^''S'o)) 

1^.(2)1 < M„ ■ ■ ■ 



where 



MJ Uist(r,ifo) 

< Af„exp{2nA(l-wc\s„(z,5o))}, z G £>, 
A := d - r + log(d(i:')/dist(r, Kq)). 



(3.29) 



Note that A is necessarily positive otherwise 6„ would be constant in D by the 
maximum principle, which is absurd. Moreover, by the regularity of 5o, it is 
known (^, Thm. 4.3.4]) that for any x e So 



lim cj£)(z. So) = 1 
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uniformly with respect to x € Sq. Thus, for any S > there exists r{6) < 
dist(5o,r) such that for z satisfying dist(z, 5'o) < r{6) we have 

l-LOD\So{^,So)<d/A. 

This, together with ( 3.29[ ), implies that for fixed S, to be adjusted later, we have 

\bn{z)\<M,,e^"', \z-xn\<r{6). 
Note that &„ is analytic in D, which, in particular, yields 



6;(^) = 



27ri 



\i-x^\=r{S) 






■d^, \z- Xn\ < r{S). 



Thus, for any z such that \z — x„| < r{5)/2 we get 
K{z)\<^-^^^^;^-^-2^r{5) 



1 4Af„e2"'5 „ _, 4Af„e2"'5 



27r r^{5) 
Now, for any x such that 
r{5) 



t{5) 



\X - Xn\ < 



p2n(5 



(3.30) 



the mean value theorem yields 

\bn{x) - bn{Xn)\ < '}j-\ N " ^n\ < -TT ■ 

r(d) 2 



Thus, for X satisfying ( 3.30[ ) and n G N2 we have 

\bn{x)\ > \bn{Xn)\ ~ \bnix) - 6„(a;„)| > Mn ~ 



M„ 



M„ 



and by (3.281 and the definition of 6„, 



<llA^) 



W2n(a;) 



> 1 -2„(d-r)^ neNa. 



(3.31) 



Now, Lemma 3.2 'a) guarantees that there exist a polynomial T of degree, 
say k, and a number j3 £ (0,7r/32) such that 



Arg(e^'^(*) Q,(t)T(f) 



< 



2/3, 



for all i G / such that {TQs){t) ^ 0, where Lp is as in (2.2). Moreover, for each 
n E N2, we choose T/^„ as in Lemma 3.2 'b) with e = S/3. Since all Ti^n are 



bounded on / by definition and have respective degrees at most I, which does 
not depend on n, there exists N3 C N2 such that sequence {Ti^n}nGN3 converges 
uniformly to some polynomial Ti on /. In particular, we have that deg(r;) < I 
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and the argument of Tijv^n lies in (—J, S) for n £ N3 large enough. Denote by 
2a the smallest even integer strictly greater than / + fc + s. As soon as n is large 
enough, since j/o G supp(j^), there exist /3i^„, . . . , f32a,n, zeros of qn,i, lying in 

{z e C : dist (z, [yo - P,yo + p]) < p} , 
such that 



2q 



EAr^ 



i=i 



1 



'^ t-^j,n 



2a 

^Ar, 



^..n) 



</3, xeM\[zjo-2p,yo + 2p]. 



Define for n e N3 sufficiently large 



_ qn{z)T{z)Ti{z) 

ll, = l(^-P,,„) 



Then 



Arg 



({PlQsqn){x)e-^^^^ 






Vln{x) 



A'8(i-wpjn_ir4jSI™-'«='*' 



< 7r/2 - (5, 



for all a; € / \ [j/o — 2p, j/o + 2/?] except if T{x)Qs{x) = 0, where 5 is chosen so 
small that 5 < (3/2. This means that for such x 



Re 



\ W2™(a;) 



> sin 5 
= sin (5 



(p;:g,(7„)(x)e^^(-) 



W2n(a;) 
92(x)Q,(a;)T(a;)Tz(a;) 



"2n(a;)n,"i(a;-/3j>) 



(3.32) 



Now, denote by rrin^i the number of zeros of qn,2 (defined in (3.18 1) of modulus 



at least 2max2:g5 \x\, and put a„ for the inverse of their product. Let to„^2 := 
deg(g„^2) — "T-n.i- Then for a; G S* we have 



(dist(5,r))"-^(l/2)"-i < |a„g„.2(x)| < ( SmaxliM (3/2)"-S (3.33) 

and since m„^i + to„^2 = deg(q„.2) is uniformly bounded with n, so is {|q;„(7„.2|} 
from above and below on S. 



Finally, ifx e S\[yo — 2p,yQ + 2p] satisfies (3.301, then by (3.311 the quantity 



in (3.321 is bounded below by 



\T{x)Qs{x)\ 

Cl 



sin(5 min^g/ |r;(x)| min^^e/ |g«,2(a;) 



^— 2nc/+2nr 



2(diam(S') + 2p)2" 
|T(a;)Q,(x)|e-2"'^+2"^ 
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where 



_ sinS mina;g/ |T;(g;)| min^g/ |Q;„gn,2(a:)P 
""' ■" 2(diam(5)+2p)2" 



by construction of T/ and (|3.33|. Thus, 
Re 



S\[vo-2p,yo+2p] 



> sin (5 



S\ly„-2p,ya+2p] 



an\'P:{t)Qs{t)qn{t) 



alP:{t)QS)qn{t) 






d\X\{t) 



V2n{t) 



d\xm 



> C\e 



-2nd+2nT 



\T{t)Q,{t)\d\\\{t) > c^e-^nd+2nir~L5) ^ (3^34) 



5n/„ 



where /„ is the interval defined by (3.301. The last inequality is true by the 



following argument. Recall that a;„, the middle point of /„, belongs to 5'o, 
where dist(S'o,2') > and Z is a finite system of points that we choose now 
to be the zeros of TQg on /, if any. Then TQg, which is independent of ti, is 



uniformly bounded below on /„ for all n large enough and (3.341 follows from 



this, the second requirement in the definition of BVT, and the fact that /„ and 
[yo — 2p, 2/0 + 2p] are disjoint for all n large enough. The latter is immediate if 
(3.261 and (3.31) are c ompar ed. 



that 



On the other hand, (3.261, applied with u„ = P^{z) /qn^2{z), and (3.33) yield 



pP„*(t)g,(t)g„(i)^d|A|(i) 



< 



cae 



-2nd 



(3.35) 



l[yo-2p,y„+2p\ V2n{t) 

This completes the proof, since S can be taken such that r — L5 > and this 
would contradict orthogonality relations (|2.15[) because, for n large enough, the 



integral in (3.34) is much bigger than in (|3.35 



Proof of Theorem\2.1\ This theorem follows immediately from Theorem|2.5|and 



the considerations leading to (2.15) 



Before we prove Theorem |2.3| we shall need one auxiliary lemma. 



Lemma 3.4. Let D be a domain in C with non-polar boundary, K' be a compact 
set in D, and {un} be a sequence of subharmonic functions in D such that 



Un{z) < M 



zeD, 



for some constant M and a sequence {£«} of positive numbers decaying to zero. 
Further, assume that there exist a compact set K' and positive constants e' and 
5' , independent of n, for which holds 



Un{z) < M 



eKnd K', cap{Kn) > S'. 
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Then for any compact set K C D\K' there exists a positive constant e{K) such 
that 

Un{z) <M -e{K), zeK, 

for all n large enough. 

Proof: Let Wn be the harmonic measure for D„ := D\ K^- Then the two- 
constant theorem f^S*, Thm. 4.3.7] yields that 

un{z) < (M-e')w„(z,i^„) + (M-e„)(l-Lj„(z,X„)) 
< Af-(e'-e„)w„(z,K„), z e D,,. 

Thus, we need to show that for any K C D\K' there exists a constant 6{K) > 
such that 

LUn{z,Kn)>S{K), Z^K. 

Assume to the contrary that there exists a sequence of points {zn}nerii C K, 
Ni C N, such that 



a'„(z„, Kn) -^0 as n ^ cx), n e Ni. 



(3.36) 



By |29l Theorem 4.3.4], LUn{-,Kn) is the unique bounded harmonic function in 
Dn such that 



hm Lu{z, Kn 



1k„(C) 

for any regular C G dDn, where 1k^ is the characteristic function of dKn- Then 



it follows from (A. 63 1 of the appendix that 



(3.37) 



where fJ-^K^.do) is the Green equilibrium measure on K relative to D. Since 
all the measures H{K„,dD) are supported in the compact set K' , there exists a 
probability measure fj, such that 

l^(K„,dD) — » /i as n — > oo, n e N2 C Ni. 

Without loss of generality we may suppose that z^ —>■ z* Cz K a,s n ~f 00, 
n G N2. Let, as usual, gD{-,t) be the Green function for D with pole a,t t E D. 
Then, by the uniform equicontinuity of {gDi',t)}teK' on K, we get 



U 



^^""■""^{zn) -^ U^{z*) 7^ as 71 ^ 00, n e N2 



Therefore, (3.36) and (3.371 necessarily mean that 



cap{Kn,dD) -^ as n ^ 00, n eN2- 



(3.38) 



By definition, l/cap(_R'„, 9D) is the minimum among Green energies of proba- 
bility measures supported on Kn- Thus, the sequence of Green energies of the 
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logarithmic equilibrium measures on Kn, fJ-K,^, diverges to infinity by (3.38). 
Moreover, since 

{g(-,i) + log|--t|},g^, 

is a family of harmonic functions in D whose moduli are uniformly bounded 
above on K' , the logarithmic energies of iik„ diverge to infinity. In other words, 

ca,p{Kn) ^ as n ^- oo, n e N2, 

which is impossible by the initial assumptions. This proves the lemma. ■ 



Proof of Theorem\2.3\ Exactly as in the proof of Theorem 2.5 we can suppose 



that all the interpolation points are contained in some compact set Kq disjoint 
from S. By virtue of the Hermite interpolation formula ( cf. [33, Lemma 6.1.2, 
(1.23)]), the error e„ := F — n„ has the following representation 



en{z) 



V2niz) 



{pn~sQsqn){t) d\{t) 



zeC\S, 



(3.39) 



{Pn~sQsqn){z) J V2n{t) Z-t' 

where Pn-s is an arbitrary polynomials in Vn-s- Since almost all of the zeros of 
qn approach S by Corollary 
in such a manner that the a 



3.3 we always can fix s of them, say ^i,n, ■ ■ ■ ,S,s 



jsolute value of 



J = l 

is uniformly bounded above and below on any given compact subset K d C\S 
for all n large enough (depending on K). In what follows we choose p„_s(z) := 
qn{z)/ls,n{z). Set also An to be 



An{z) 



i\pl^s\Qsls,nm dX{t) 



V2n(t) Z-t 

First, we show that 

|A„|i/2"^APexp{-c((7,C\5)} 



, zc£\S. 



(3.40) 



on compact subsets of C\5, where c((7, C\S') is defined in (A.6OI of the appendix 



Clearly, for any compact set K <z€\S there exists a constant c{K) , independent 
of n, such that 



\An{z)\<c{K) 



pI- 



W2ri 



z e K. 



(3.41) 



Let Vn and cr„ be the counting measures of zeros of p^j_^ and U2n, respectively. 
Then 



lim sup 



V2n{t) 



l/2r, 



= lim infexp {[/'"" "''"(i)} = exp iu'^-" {t)\ 
= exp{-c(cr;C\S')} q.e. on S (3.42) 
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by Theorem 2.1 the lower envelope theorem [3DJ Thm. 1.6.9], and (A. 60 1 of the 
appendix. Moreover, by the principle of descent fSO', Thm. 1.6.8], we get that 



lim sup 



pl-At) 



l/2r: 



< exp 



(^-^{1)1 = exp{-c(CT; C \ S)} (3.43) 



V2n(t) 

uniformly on S, where the last equality holds by the regularity of S. Now it is 



immediate from (3.421 and (3.43) that 



lim 

n — >oo 



Pn — s 



V2n 



l/2n 



= exp{-c(cr;C\S')}. 



(3.44) 



Indeed, since the whole sequence {vn} converges to ct, (3.421 holds for any 
subsequence of N. Thus, there cannot exist a subsequence of natural numbers 



for which the limit in (3.441 would not hold. Suppose now that (3.401 is false 



Then there would exist a compact set K' C C\ S and e' > such that 

capLeK': |A„(z)|i/2"-exp{-c(o-;C\S')} >e'}74 0. (3.45) 



Combining (3.45), (3.44), and (3.41) we see that there would exist a sequence 



of compact sets Kn C K' ,ca.p{Kn) > 5' > 0, such that 
|A.(^)r/'" < exp{~c{a; C\S)}~ e', z e X„. 



(3.46) 



Now, let r be a closed Jordan curve that separates S from Kq, a compact set 
containing all the zeros of V2n, and K'. Assume further that S belongs to the 
bounded component of the complement of T. Observe that (l/2n) log |A„| is a 
subharmonic function in C \ S*. Then by ( |3.41 1, (3.44), and (3.46) enable us to 
apply Lemma 3.4 with AI — ~c{a;C\S) which yields that there exists e(r) > 
such that 



|^n(^)r/'" < exp{-c(a; C \ 5) - e(r)} 
uniformly on T and for all n large enough. Define 

dz 



(3.47) 



^n 



Ti{z)T{z)U^^{z)A^{z) ^ . 

p ZTTJ 



where the polynomials T; and T are chosen as in Theorem 2.1 (see discussion 
after ( |3.31[ )). We get from ( |3.47[ ) that 

(3.48) 



limsup jy^" < exp{-c(a;C\ S*) - e(r)}. 



In another connection, Fubini-ToncUi theorem and the Cauchy integral formula 
yield 

„2 



^n 



^ ' ^ ' '^ ' \J V2n{t) z-t) 2m 

\ql{t)\^ATQs){t)e^'^^'U\X\{t) 

V2n{tj 



(3.49) 
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Exactly as in (3.321, we can write 



Re — > sm{S) 



{TiTQsm 



te L 



(3.50) 



V2n{t) 

where / is the convex hull of S and (5 > has the same meaning as in the proof 



of Theorem 2.1 (see construction after (3.29l). Thus, we derive from (3.49) and 



( |3.50| that 

J„ > sin(,5) / |6„(t)| \{TiTQ,){t)\d\X\{t), &„ := ql/v2n- 



(3.51) 



Let 5*0 be a closed subset of S of positive capacity that lies at positive distance 



from the zeros of TQg on / (as in the proof of Theorem 2.1 we refer to [2J [3] 
for the existence of this set). Further, let a;„ G Sq be such that 

ll^nllso = \bn{Xn)\- 



Then it follows from (3.421 that 



\bn\\s„ > exp{-2n(c(CT; C\S) + e)} 



for any e > and all n large enough. Proceeding as in Theorem 2.1 (see 



equations (3.301 and (3.31)), we get that 



\bn{t)\ > - exp{-2n(c(a; C\S) + e)}, t £ /„, 



(3.52) 



where 



/„:={a;e5o: |x-a;„| < r^e-^"^} 

and rg is some function of 5 continuous and vanishing at zero. Then by combin- 
ing (3.51 ) and (3.52), we obtain exactly as in (3.34) that there exists a constant 



ci independent of n such that 



J„>sin(5) / |6„(t)| \{TiTQ,){t)\d\X\{t)>cie^^{-2n{c{a-C\S) + e + L5)}. 

J Ir, 

Thus, we have that 

liminf jy2« > exp{-c(CT; C \ 5) - e - L5}. (3.53) 

n — >oo 

Now, by choosing e and 5 small enough that e + L5 < e(r), we arrive at contra- 



diction between (3.48) and (3.53). Therefore, the convergence in (3.40) holds 



Second, we show that 



V2n{z)ls,n{z) 



ql{z)Qs{z) 



l/2n 



3'exp{c(a;C\5)-C/^\5(z)} 



(3.54) 
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on compact subsets of C \ S*. Let K C C\S he compact and let [/ be a bounded 
open set containing K and not intersecting S. Define 

QnA^) '■= H (z ~ and QnA^) -^ qn{z)/qn.l{z). 

ieu: g,.(C)=0 



Corollary 3.3 yields that there exists fixed m G N such that deg(g„_i) < m. 
Then 



V2n{z) 



Tnai^) 



l/2n 



exp 



{c/^--(z)}=exp{c(a;C\5)-C/^\5(z)} 



uniformly on K by Theorem 2.1 definition of U2n, and (A. 62 1 of the appendix. 



Moreover, it is an immediate consequence the choice of ls,m the uniform bound- 
edness of the degrees oi q^^Qs, and ^29j Thm. 5.2.5] {ca.\i{{z : \{qn,iQs){z)\ < 

e}) = edcg(gS,iQ,)) ^^^^ 

l/2n 



^s,n \Z ) 



Qliiz)QM 



z€K. 



and (3.54) 



Thus, we obtain (3.541. It is clear now that (2.9) follows from (3.39), (3.40) 



Proof of Theorem 2.4' Inequality (2.11) is trivial for any rj £ S' r\ S. Suppose 



now that rj £ S' \ S and that 771(77) < "^('y)- This would mean that there exists 
an open set U, U H S — {?/}, such that 777(77, U) < m{ri) and therefore would 
exist a subsequence Ni C N such that 

#{5„nc/}< 777(77), neni. 



It was proved in Theorem |2.3| that {n„} converges in capacity on compact sub- 
sets of C \ 5" to -F. Thus, {n„}„gNj is a sequence of meromorphic (in fact, 
rational) functions in U with at most 777(77) poles there, which converges in ca- 
pacity on C/ to a meromorphic function F\jj with exactly one pole of multiplicity 
m(r]). Then by Gonchar's lemma jT4} Lemma 1] each n„ has exactly 777(77) poles 
in U and these poles converge to 77. This finishes the proof of (2.11 1. 



3.3 



Now, for any r] € S' \ S the upper characteristic 7(7(77) is finite by Corollary 



Therefore there exist domains D^i, D^ n S 



{77}, such that 777(77) 



rn{r],Dj^), 77 G 5" \ S. Further, let 9{-) be the angle function defined in (2.10) 
for a system of 777 intervals covering S and let S'„ = {^i.n, • • ■ , ^d„,n}- Then by 
Lemma [3.11 we have 



^{n-e{i,^n)) + {n-dn)TT<V{ip) + VA + {m+s-l)7T+y^m{C)e{0. (3.55) 



i=i 



Ces' 



Then for 77 large enough (3.55) yields 



E E (^-^(a"))-m(77)(^-%)) <y, 



ves'\s Uj,„G-D, 
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where V was defined in (2.13). Thus, 



< Enes'\s *{Sn n D,} (max^eA, ^(0 - 0{v)) +V 

for all n large enough. However, since {max„>Ar ^{S„r\D,-i}}N^f!i is a decreasing 
sequence of integers, 771(77) — ^('7? ^r;) = #{'5'nnl?,,} for infinitely many 77 £ N. 



Therefore, we get from (3.561 that 



J2 iHv) - Hv)) i^-div)) <V+ Y. ^(^) ( S|f ^(^) - ^(^) ) ■ (3-57) 
Observe now that the left-hand side and the first summand on the right-hand 



side of (3.571 are simply constants. Moreover, the second summand on the 



borhoods D„. Thus, (2.12) follows 



right-hand side of (3.57) can be maid arbitrarily small by taking smaller neigh- 



4 Numerical Experiments 

We restricted ourselves to the case of classical Pade approximants and we con- 



structed their denominators by solving the orthogonality relations (2.15) with 
W2n = 1- Thus, finding these denominators amounts to solving a system of lin- 
ear equations whose coefficients are obtained from the moments of the measure 
A. 

In the numerical experiments below we approximate function F given by the 
formula 

F{z) = 7 ---(3 + 7) ' / 



- , 6/7,-1/8] ^-'t 7(2/5,1/2] < - 27 Z-t 

+(2-47)/ '^"^ 

i[2/3,7/8] ^ ^ * 

1 2 6 



(z + 3/7-4i/7)2 (z - 5/9 - 3i/4)3 {z + l/b + Qi/lY' 

On Figures la and 2a the solid lines stand for the support of the measure, 
diamonds depict the polar singularities of F, and disks denote the poles of the 
corresponding approximants. Note that the poles of F seem to attract the 
singularities first. On Figures lb and 2b the absolute value of the error on 
the unit circle is displayed for the corresponding approximants. The horizontal 
parts of the curves are of magnitude about 10^'^ on Figure lb and of magnitude 
about 10^^ on Figure 2b. 

5 Appendix 

Below we sketch some basic notions of logarithmic potential theory that were 
used throughout the paper. We refer the reader to the monographs [29l [30] for 
a complete treatment. 
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0.2 0.2 0.4 0.6 



Figure 1: Poles of Tin and the error \F — Ilui] on T 



The logarithmic potential and the logarithmic energy of a finite positive mea- 
sure fx, compactly supported in C, are defined by 



t/^(z):= / log j-^dfi{t), ze 



and 



/[m] := / t/^(z)dM(z) - / / log-— dfi{t)dfi{z) 



(A.58) 



(A.59) 
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Figure 2: Poles of 1120 the error |F — n2()| on T 



respectively. The function V^ is superharnionic with values in (— oo,+oo], and 
is not identically +oo. It is bounded below on supp(/i) so that I[fj] G (— oo, +oo]. 

Let now _E C C be compact and A{E) denote the set of all probability 
measures supported on E. If the logarithmic energy of every measure in A(E) 
is infinite, we say that E is polar. Otherwise, there exists a unique iie € A(E) 
that minimizes the logarithmic energy over all measures in A{E). This measure 
is called the equilibrium distribution on E. The logarithmic capacity, or simply 
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the capacity, of E is defined as 

cap(£;) = exp{-/[^B]}. 

By definition, the capacity of an arbitrary subset of C is the supremum of the 
capacities of its compact subsets. We agree that the capacity of a polar set is 
zero. We define convergence in capacity as follows. We say that a sequence of 
functions {ft.„} converges in capacity to a function h on a, compact set K if for 
any e > holds 

ca,p{{z£K: \{hn — h){z)\ > e}) ~^ as n — > oo. 

We also say that a sequence converges in capacity in an open set f2 if it is 
converges in capacity on any compact subset of 51. 

Another important concept is the regularity of a compact set. We restrict 
to the case when E has connected complement, say il. Then E is called regular 
if the Dirichlet problem on dH. is solvable, in other words, if any continuous 
function on dfl is the trace (limiting boundary values) of some function harmonic 
in fi. Thus, regularity is a property of dfl rather than E itself. It is also known 
[501 pg. 54] that E is regular if and only if f/'^^ is continuouaj in C 

Often we use the concept of balayage of a measure ([30, Sec. 11.4]). Let D be 
a domain (connected open set) with compact boundary dD whose complement 
has positive capacity, and /i be a finite Borel measure with compact support in 
D. Then there exists a unique Borel measure fi supported on dD, with total 
mass is equal to that of /i, whose potential U^ is bounded on dD and satisfies 
for some constant c{fi; D) 

U^iz)^U^'iz) + c{fi;D) forq.e. zeC\D. (A.60) 

Necessarily then, we have that c{fj,; D) = if D is bounded and c(/j,; D) = 
J gD{t,oo)dfj,{t) otherwise, where f/£i(',oo) is the Green function for D with 



pole at infinity. Equality in (A.60 1 holds for all z S C \ D and also at all regular 
points of dD. The measure /I is called the balayage of /i onto dD.It has the 
property that 

hdfi^ j hdjl (A.61) 

for any function h which is harmonic in D and continuous in D (including at 
infinity if D is unbounded). From its defining properties 'jl has finite energy, 
therefore it cannot charge polar sets. 

In analogy to the logarithmic case, one can define the Green potential of a 
positive measure ^ supported in a domain D with compact non-polar boundary. 



The only difference is now that, in (A. 58 1, the logarithmic kernel log(l/|z — t\) 



gets replaced by gD{z,t), the Green function for D with pole at t € D. The 



^Since supp(/i£) C dO, |30l Cor. 1.4.5], it is again enough to check continuity of U^^ only 
on dVl. 
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Green potential relative to the domain D of a finite positive Borel measure /i 
compactly supported in D is given by 

U^{z)^ J gD{z,t)dfi{t). 

It can be re-expressed in terms of the logarithmic potentials of fi and of its 
balayage Jl onto dD by the following formula [301 Thm. II. 4. 7 and Thm. II. 5.1]: 



[/?-A'(^) = dfi; D) - C/^(z), z£D, (A.62) 



where c{fj,;D) was defined after equation (A. 60 1. Moreover, (A.62 1 continues to 
hold at every regular point of dD; in particular, it holds q.e. on dD. 

Exactly as in the logarithmic case, if E' is a compact nonpolar subset of D, 
there exists a unique measure ^i(E.dD) G A(-^) that minimizes the Green energy 
among all measures in A(ii^). This measure is called the Green equilibrium 
distribution on E relative to D. In addition, the Green equilibrium distribution 
satisfies 

t/r-W-^;^^^^, forq.e. zeE, (A.63) 

where cap(ii', dD) is Green (condenser) capacity of E relative to D which is the 
reciprocal of the minimal Green energy among all measures in A(£'). Moreover, 
equality in (A.63[) holds at all regular points of E. 
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